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Abstract. The annual cycle of the net mass transport across the extratropical tropopause 
is examined. Contributions from both the global-scale meridional circulation and the mass 
variation of the lowermost stratosphere are included. For the northern hemisphere the mass 
of the lowermost stratosphere has a distinct annual cycle, whereas for the southern hemisphere, 
the corresponding variation is weak. The net mass transport across the tropopause in the 
northern hemisphere has a maximum in late spring and a distinct minimum in autumn. This 
variation and its magnitude compare well with older estimates based on representative 90 Sr 
mixing ratios. For the southern hemisphere the seasonal cycle of the net mass transport is 
weaker and follows roughly the annual variation of the net mass flux across a nearby isentropic 
surface. 


1. Introduction 

Soon after the start of nuclear bomb testing in the strato- 
sphere in the late 1950s it was found that in the northern 
hemisphere the near surface value of radioactive debris showed 
a spring maximum and an autumn minimum fe.g., Fry et ai, 
I960]. It was argued [e.g., Staley , 1962] that this seasonal 
variation in radioactivity could be linked to a seasonal cycle 
in stratosphere-troposphere exchange (STE). Detailed case 
studies supported this idea [e.g., Danielsen , 1968] and indi- 
cated that upper level baroclinic events, such as tropopause 
folds, are important for mass exchange. Early studies of the 
seasonal variation of near-tropopause cyclogenesis [ Mahlman , 
1969] and more recent studies of the frequency of cutoff lows 
[Price and Vaughan , 1992] and tropopause folds [Beekmann 
et al, 1996] however showed no preferred occurrence in the 
spring season. This suggests that other factors than the number 
of events must be important for producing the apparent spring 
maximum and autumn minimum in exchange. 

In a recent review, Holton et al. [1995] emphasized the global- 
scale aspects of exchange. They argued that wave-induced 
zonal forces in the extratropical stratosphere induce a global- 
scale meridional circulation in which mass is pulled upward and 
poleward in the tropics and pushed downward in the 
extratropics. This global-scale meridional circulation, which 
must be accompanied by radiative heating and cooling, is 
often referred to as the diabatic circulation [ Brewer , 1949; 
Dobson , 1956]; it shows a clear seasonal cycle in both 
hemispheres [Rosenlof and Holton , 1993; Rosenlof 1995]. 
The maximum downward mass flux across a control surface that 
lies in the lower stratosphere (e.g., at 100 or 70 hPa) occurs, 
however, during the winter and not in late spring. 

To understand mass and tracer transport into and out of the 
middle stratosphere the replacement of the tropopause by a 
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more convenient isobaric or isentropic control surface located 
in the lower stratosphere is useful. However, for some pur- 
poses the mass transport across the actual tropopause is 
required. The net downward mass fluxes across the extra- 
tropical tropopause and across an isobaric or isentropic 
control surface that nearly coincides with the tropopause in the 
tropics should be equal for a sufficiently long time average 
[Holton et al. , 1995]. Such equality will not hold on seasonal 
or shorter timescales since the amount of mass in the layer 
between the extratropical tropopause and an isobaric or 
isentropic surface coinciding with the tropical tropopause may 
vary with time. 

It has long been recognized that stratospheric mass varia- 
tions due to seasonal tropopause height variations can contri- 
bute to stratospheric tropospheric air exchange [Staley y 1962; 
Reiter , 1975; Robinson , 1980; Murgatroyd and O'Neill , 1980]. 
Because of the lack of global tropopause analyses the global- 
scale effects of these variations have only been roughly 
estimated. In this paper the contributions of the seasonal 
variations of the mass of the lowermost stratosphere and of the 
global-scale meridional circulation in the stratosphere to the 
variations in net mass transport across the extratropical 
tropopause will be explored. 

2. Model for Net Mass Transport 

The simple model used is shown in Figure 1. The shaded 
region corresponds to the “lowermost stratosphere,” also 
referred to as the stratospheric middle world [ Hoskins , 1991; 
Holton et al. , 1995]. In this region mass can be exchanged 
with the troposphere along isentropic surfaces (dashed-dotted 
curves). It is separated by a potential vorticity surface (bottom 
solid curve) from the troposphere and by an isentropic surface 
(top solid curve) from the overworld, the region of the strato- 
sphere where mass can only be exchanged by crossing isen- 
tropic surfaces. The global-scale meridional circulation (open 
arrows) supports a mass flow (F in ) directed into the lowermost 
stratosphere across the upper boundary. By mass continuity 


15,071 



15,072 


AFPENZELLER ET AL.: CROSS-TROPOPAUSE MASS TRANSPORT 



Pole mid-latitude equator 

Figure 1 . The simple hemispheric net mass transport model. 
See text for details. 


the flow out of the lowermost stratosphere into the tropo- 
sphere (F out ) is simply the difference between the inflow from 
the overworld and the rate of increase ot mass ( M ) of the 
lowermost stratosphere: 

E out (0 - ~M + VO (1) 

dt 

where F in and F oul are both defined to be negative for downward 
directed mass fluxes. 

In this model any two-way transport ( F between the tropo- 
phere and the lowermost stratosphere that is characterized by a 
zero net mass transport across the tropopause is not included. 
Assuming that M and Fj n can be calculated from “observed” 
data, equation (1) provides a simple way to estimate the net 
mass flux F out across the tropopause, without doing any flux 
calculation at the tropopause itself. 


3. Mass of Lowermost stratosphere 

The mass of the lowermost stratosphere is given by 



1 

— dp cos 0 dtp dK 
8 


(2) 


where g is gravity, A is longitude, <p is latitude, p is pressure, 
and P\ (A, <p) and F 2 (A, <j>) represent the pressures on the lower 
and upper boundary surfaces, respectively. The former is taken 
to be the 2 PVU potential vorticity surface (1 PVU = 10 -6 K m 2 
s -1 kg -1 ), while the latter is taken to be the 8 - 380 K 
isentrope [ Hoskins , 1991; Reid and Gage , 1981; Hoskins et 
al ., 1985; Holton, et al, 1995]. The global tropopause is 
therefore defined by a combination of a potential vorticity and 
a potential temperature surface. The sensitivity of the calcula- 
tion to the chosen PV and 0 value is discussed below. 

The atmosphere’s PV and 8 structure was calculated for two 
calendar years 1992 and 1993 from daily (00:00 UT) global 
stratospheric analyses prepared by the United Kingdom 
Meteorological Office (UKMO) for the Upper Atmosphere 
Research Satellite (UARS) project [Swinbank and O'Neill , 
1994]. Horizontal winds, geopotential heights, and tempera- 
ture are provided with a horizontal resolution of 2.5° latitude 
by 3.75° longitude, and vertical resolution ~50 hPa in the 
lowermost stratosphere. 

According to equation (2) any variation of the mass of the 
lowermost stratosphere is linked to a net variation of the 
pressure distribution on the tropopause or the 380 K surface. 
An example of the instantaneous pressure distribution for both 
surfaces is given in Figure 2. For both charts the same contour 
interval of 25 hPa is used. Pressure variations on the isen- 
tropic surface are relatively small and nearly zonally sym- 
metric; pressure variations on the tropopause are large and 
have strong deviations from zonal symmetry. The latter are 
linked to baroclinic developments and subsynoptic scale 



Figure 2. Instantaneous view of the pressure distribution on the northern hemisphere tropopause (a) and 
380 K isentrope (b). Contour interval is 25 hPa for both panels, day number 248 of year 1993. 
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Figure 3. Zonal mean position of the upper and lower 
boundaries of the lowermost stratosphere for January (solid) 
and July (dashed) and annual mean (thin long dash). Data are 
averaged in pressure coordinates and the monthly mean is 
based on zonal mean daily data for the year 1993. 

disturbances, such as cutoff systems and stratospheric intru- 
ions [e.g., Hoskins et ai, 1985; Appenzeller et ai , 1996]. 

The monthly mean positions of the tropopause and the 
380 K isentrope based on zonal mean daily data are shown in 
Figure 3 for January and July 1993. In the tropical region the 
pressure of the 380 K isentrope appears to be relatively con- 
stant with a typical value around 100 hPa and an annual varia- 
tion of less than 10 hPa (Figure 4a). In the extratropical 


regions the mean height of the 380 K surface is for both hemi- 
spheres at a maximum in winter and at a minimum in the 
summer season. This up and down oscillation is consistent 
with the seasonal variation in the strength of the stratospheric 
polar vortex. The relative amplitude (Figure 4a) of the 
oscillation is slightly stronger for the southern hemisphere 
than for the northern hemisphere, but the overall pattern is 
comparable. 

In the northern hemisphere the tropopause moves upward 
and poleward during the transition from winter to summer 
(Figure 3, Figure 4b). In the southern hemisphere a similar 
shift in the mean position of the tropopause occurs only 
between 30°S and 40°S. In more poleward southern latitudes 
the tropopause moves downward and equatorward during the 
transition from winter to summer. Thus the seasonal cycle in 
the southern hemisphere tropopause pressure is a “seesaw” 
oscillation in which a relatively high tropopause in southern 
subtropical regions is accompanied by a relatively low tropo- 
pause in poleward latitudes and vice versa (Figure 4b). Note 
that at any given time or location the position of the tropo- 
pause may differ substantially from the seasonal meridional 
mean. To incorporate this variability and to avoid any bias 
due to the averaging procedure, the mass variation of the 
stratospheric middle world was calculated from equation (2) 
using daily global data. 

The calculated mass of the lowermost stratosphere has a 
distinct seasonal cycle (Figure 5a). For the northern hemi- 
sphere a maximum occurs in midwinter and a minimum in late 
summer, whereas the corresponding seasonal variation for the 
southern hemisphere is much weaker in amplitude. Some 
higher frequency variations were also found. For the northern 
hemisphere their amplitudes are an order of magnitude smaller 
than the amplitude of the annual variation. 

The mass of the stratospheric overworld (Figure 5b) also 
shows a seasonal cycle. In this case the amplitudes of the 
cycles in the two hemispheres are comparable, and the phases 



Figure 4. Annual variation of the zonal mean relative pressure on the 380 K isentrope (a) and the 
tropopause (b) for the two years 1992 and 1993. The pressure is given relative to the annual zonal mean (as 
defined in Figure 3) with a contour interval of (a) 5 hPa and (b) 10 hPa, negative values dashed, zero line is 
dashed-dotted. 
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Figure 5. Annual variation of the mass (I0 ,7 kg) for 
various subregions of the stratosphere for the two years 1992 
and 1993. Plain lines for the northern hemisphere (NH) and 
dashed lines for the southern hemisphere (SH). (a) Lowermost 
stratosphere, (b) Stratospheric overworld, (c) Total strato- 
sphere. 


are opposite to those of the lowermost stratosphere, consis- 
tent with the observed up and down oscillation of the 380 K 
surface (Figure 3 and 4a). The global stratospheric overworld 
experienced a weak mass loss of -1.5% during the calendar 
year 1992 and 1993, as was also noted in a recent analysis by 
Yang and Tung [1996]. 

Combining Figures 5a with 5b results in the mass variation 
of the total stratosphere (Figure 5c). For the northern hemi- 
sphere this closely follows the variation of the mass of 
the lowermost stratosphere, which implies that the latter is 
mainly linked to the seasonal variation of the mean position 
of the tropopause. For the southern hemisphere the seasonal 
cycle of the total stratospheric mass is very weak. This is 
consistent with the seesaw nature of the tropopause pressure 
oscillations and suggests that the stratospheric air mass is 
shifted from low latitudes to more poleward latitudes and vice 
versa without changing the hemispheric averaged tropopause 
height. The pronounced asymmetry between the two hemi- 
spheres that we have found in the annual cycle in total 
stratospheric mass stands in contradiction to earlier studies 
that attributed the seasonal mass variation to a mass shift 
between the two hemispheres [ Reiter , 1975] or suggested only 
a slightly weaker seasonal variation for the southern hemi- 
sphere [Murgatroyd and O'Neill , 1980]. 

The montly mean contribution to F out due to mass changes 
of the lowermost stratosphere as given by dM/dt in equation 
(1) was calculated on the basis of monthly averages of the 
daily mass values. For the northern hemisphere a maximum 
downward flux (negative dM/dt) occurred in late spring 
(Figure 6), whereas the seasonal cycle for the southern hemi- 
sphere is less distinct. Also shown are similar calculations 
using different definitions for the tropopause and the upper 



time [month] 

Figure 6. Annual variation of the net mass flux across the 
tropopause due to the mass variation of the lowermost strato- 
sphere alone, using various definitions of the boundaries for 
the two years 1992 and 1993, Solid curve, 2 PVU and 380 K; 
dotted curve, 3 PVU and 380 K; and dashed curve, 2 PVU and 
370 K. (Values are given in 10 9 kg s _J ; negative values denote 
mass flux out of lowermost stratosphere.) 
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boundary of the lowermost stratosphere. Although the actual 
total mass values differ, the mass fluxes show only a very weak 
sensitivity to the choice of control surfaces. This indicates 
that the above described seasonal variations and their asym- 
metry between the two hemispheres are quite robust and should 
also be found using higher-resolution data. 

In the northern hemisphere the annual mass variation is 
- 0.6 x 10 17 kg for the total stratosphere (i.e., -10 % of the 
mass of the stratosphere) and - 0.75 x 10 17 kg for the lower- 
most stratosphere (i.e., -30 % of the mass of the lowermost 
stratosphere). The mass variation of the total stratosphere can 
be used to estimate the minimal annual mass transport across 
the tropopause. The value given above is of the same order of 
magnitude as suggested in earlier studies and is about 60% of 
that found by Robinson [1980], who used a crude estimate 
based on the jet displacement alone. 


4. Mass Transport Across Upper Boundary 

To calculate the mass flux Fj n across the upper boundary of 
the lowermost stratosphere, the mass flux across the 380 K 
isentrope due to the global-scale meridional circulation was 
estimated. In isentropic coordinates the vertical velocity is 
directly given by the net diabatic heating rate Q. The zonal 
mean mass flux rate can then be approximated by [Tung, 1982; 
Pawson and Harwood , 1989] 


with 


F° ~ cos (<p) a Q 


(3) 


a - 



(4) 


representing the zonal mean of the density in isentropic 
coordinates. The monthly averaged net diabatic heating rates 
Q were calculated from the radiative algorithm of Olaguer 
et ai [1992] with UARS-based constituent measurements and 
UKMO temperatures as input. The ozone measurements were 
taken from the Microwave Limb Sounder (MLS), whereas water 
vapor and methane measurements were from the Halogen 
Occultation Experiment (HALOE) , except for the year 1992 
when H2O was taken from MLS. At lower pressure levels, 
climatological values were used. The heating rates were 
calculated on pressure levels and were corrected for mass 
imbalances to assure that the net vertical mass flux across a 
pressure surface is zero [Rosenlof, 1995; Shine, 1989; 
Solomon et ai, 1986]. The mass fluxes calculated from 
equation (3) using corrected heating rates are nearly equivalent 
to the mass fluxes obtained using radiatively derived trans- 
formed Eulerian mean (TEM) velocities, such as those of 
Rosenlof [1995], when the latter are interpolated onto 
isentropic surfaces, and correction is made for the temporal 
variation of the zonal mean pressure on the isentropic surface 
(see Appendix). 

For the northern hemisphere the extratropical net downward 
mass flux (Fj n ) across the 380 K surface has a maximum in 
midwinter and a minimum in late spring to early summer 
(Figure 7a, plain thick curve). This variation is consistent 
with the idea that the annual cycle in zonal momentum forcing 
in the stratosphere controls the downward mass flux into the 
lowermost stratosphere [Rosenlof 1995; Holton et ai , 1 995]. 

The amplitude of the seasonal cycle of F in is much weaker 
when evaluated on an isentropic surface than when evaluated 


on an isobaric surface with the same annual mean net mass 
flux. In the southern hemisphere the amplitude evaluated on 
an isentropic surface is much closer to that evaluated on an 
isobaric surface (Figure 7, dashed curves). In both hemispheres 
the phase of the annual cycle in cross-isentropic flux leads that 
of the cross-isobaric flux by roughly 1 month. The corre- 
sponding tropical net upward mass flux across the 380 K sur- 



time [month] 



time [month] 

Figure 7. (a) Annual variation of the extratropical net 
downward mass transport across the 380 K isentropic surface 
(F380, thick curves; plain for NH and dashed for SH) and the 
118-hPa pressure surface (Fug, thin curves; plain for NH and 
dashed for SH) for the two years 1992 and 1993. (b) Same as 
Figure 7a but for the tropical net upward mass transport (thick 
plain across 380 K, thick dashed across 1 18 hPa). (Values are 
given in 10 9 kgs _1 ; negative values denote downward mass 
flux.) 




15,076 


APPENZELLER ET AL.: CROSS-TROPOPAUSE MASS TRANSPORT 


face is shown in Figure 7b (solid curve). It has a negligible 
seasonal cycle in comparison to the one across a pressure 
surface (dashed curve). 

The total annual mean net downward mass fluxes are sum- 
marized in Table 1. It is important to bear in mind that these 
estimates are subject to substantial uncertainties. In the lower 
stratosphere the net radiative heating rates are a small residual 
of several large terms and there exists no physically based 
method to correct the heating rates to achieve global mass 
balance [Shine, 1989]. An estimate of the uncertainty in the 
corresponding mass flux rates can be gained by comparing the 
annual mean net downward mass flux as given in Table 1 with 
the one calculated in a recent study by Yang and Tung [1996]. 
They used a slightly different version of the same radiative 
transfer code [Olaguer et ai, 1992] together with a two- 
dimensional isentropic model in assimilation mode to esti- 
mate 5-day-averaged mass fluxes across a 385 K surface. For 
the northern hemisphere their annual averaged downward mass 
flux is ~ 20% lower than the one calculated here. For the 
southern hemisphere the difference is ~ 50% (Table 1). The net 
mass flux calculated here was also sensitive to the 6 surface 
selected (40% higher at 370 K and 25% smaller at 390 K). 
However, the amplitude of the seasonal variation was 
comparable. 

5. Mass Transport Across the Tropopause 

Substituting Fj n and dM/dt in equation (1) yields a simple 
estimate of the net mass transport F oul across the tropopause. 
In the northern hemisphere this transport has a distinct 
seasonal cycle (solid curve in Figure 8a) with a maximum in 
late spring and a secondary maximum in midwinter. A distinct 
minimum occurs during autumn. Compared to the mass flux 
across a 380 K surface (thin dashed curve in Figure 8), the 
annual amplitude variation is stronger and the maximum occurs 
with a phase lag of approximately 5 months. For the southern 
hemisphere (Figure 8b) the seasonal cycle was much weaker. 
The mass transport has a maximum in midwinter and is domi- 
nated by the mass flux across the 380 K surface. 

6. Conclusion 

The global-scale meridional circulation and the mass 
variation of the stratosphere have been analyzed to estimate 
the net mass transport across the tropopause. The mass 
calculations and the tropopause analyses were based on the 
atmosphere's potential temperature and potential vorticity 


Table 1 . Total Annual Downward Mass Flux Across the 380 
and 385 K Isentropic Surface in Units of I0 17 kg/year 



380 K 

385 K 

Y&T385 K 

Northern hemisphere 

3.5 

3.0 

2.4 

Southern hemisphere 

3.3 

2.8 

1.4 


The values represent averages over two years (1992 - 1993). 
The values in column 3 are taken from Yang and Tung [1996] 
and are averaged over 3 years (1992 - 1994). 


structure, using assimilated UKMO analysis data. In the 
northern hemisphere the annual zonal mean changes of the 
tropopause height are linked to a net mass variation of the 
lowermost stratosphere. The mass flux due to this variation 
substantially contributes to the net mass flux across the 
tropopause and has a maximum in late spring and a minimum 
in early autumn. For the southern hemisphere the corre- 
sponding annual variation is weak, even though the mean 
position of the tropopause shows a seasonal variation. The 



time [month] 



time [month] 

Figure 8. Annual variation of the net downward mass trans- 
port across the extratropical tropopause (F ou{ , thick plain 
curve), the 380 K isentropic surface (F^o* thin dashed curve), 
and mass flux according to 90 Sr measurements (Fg qs m dotted 
line, see text). Shown are the two years 1992 and 1993. (a) 
Northern hemisphere, (b) Southern hemisphere. (Values are 
given in 10 9 kgs _l ; negative values denote downward mass 
flux.) 
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global scale meridional circulation was approximated by radia- 
tively derived residual velocities, which were calculated based 
on UARS constituent measurements. The maximum downward 
mass flux across the upper boundary of the lowermost 
stratosphere (the 380 K isentrope) occurs during the winter 
season, and the minimum flux occurs in early summer. This 
cycle has a similar phase but weaker amplitude than the annual 
cycle in the downward mass flux computed for a nearby 
pressure surface. 

According to our simple model shown in Figure 1 it is the 
combined effects of the seasonal variation of the global-scale 
meridional circulation and the seasonal variation of the mass 
of the lowermost stratosphere that determine the net mass 
transport across the tropopause into the troposphere. For the 
northern hemisphere this flux shows a late-spring maximum 
and an autumn minimum. In the southern hemisphere the 
corresponding annual variation is dominated by the mass flux 
across the upper boundary of the lowermost stratosphere. The 
annual variation and its magnitude compare well (Figure 8a) 
with a northern hemisphere estimate based on representative 
90 Sr mixing ratios [Danielsen and Mohnen , 1977, equation 
(I)]. The mass transport cannot be directly interpreted in terms 
of tracer fluxes; but it can be anticipated that tracer inputs into 
the troposphere will be strongly correlated with the observed 
variations in the net mass transport. 

Appendix 


Here, we used 



(A5) 


where the subscripts of the brackets^ denote the particular 
vertical coordinate, s = y and A = 6 or \ff . 

By comparing equations (A4) with (A3), it follows that 


cos (0) p 0 Q 


30 ( 

dz [_v (b )q 


COS (0) p a 



(A6) 


where the partial time derivative of 0 is also transformed 
according to equation (A5) with s = /. Noting that 
z t = ~(H/p)p t and 6, = ~(p/H)0 p where// repre- 
sents the scale height, p the pressure and the subscripts denote 
partial derivatives, equation (A6) can be rewritten as 


_ 30 

cos (<p)p 0 Q = — 
dp 


P_( dy 1 


P o COS(0)| 


[ 

dt 



■ (A7) 


Finally, approximating the zonal mean of the isentropic 
density as a ~ - (l/g)(op/<90) and using p/Pq ~P s /po 
where p s and p 0v represent the constant surface pressure and 
surface density, it follows that 


The following derivation shows that the zonal mean mass 
flux rate across an isentropic surface as given in equation (3) is 
closely related to the mass flux computed from radiatively 
derived residual velocities (i.e., transformed Eulerian-mean 
(TEM) velocities [Andrews and McIntyre , 1976; Andrews et al . , 
1987]). By using the mass streamfunction \]f defined as 
[Rosenlof and Holton , 1995; Rosenlof “ 1995] 


3 

— \y = - cos (0)p n v 
3z. 

3 

— y / = cos (0)p n w 
3y 


(Ala) 


(Alb) 


the approximate TEM thermodynamic equation 


30 _* 30 * 30 

— + v — + w — ~ Q 
3t 3y 3z 


can be rearranged to 

3\f/ 30 3y/ 30 

= cos (0) p 0 

3y 3z 3z 3y 



30 ' 
3t _ 


(A2) 


(A3) 


p s (&¥) 1 (dp\ 

cos (0) a Q = 1 — — + cos (0) I — I . (A8) 

//p 0 g V 3y )q g^3tJo 

Here, it is also assumed that the zonal mean net diabatic 
heating rate in isentropic coordinates can be approximated by 
Q averaged in log-pressure coordinates. 

The zonal mean mass flux across an isentropic surface is 
therefore given by two contributions; the first term in equation 
(A8) represents the mass flux due to the residual velocity 
across the isentropic surface fixed in space, whereas the second 
term represents the contribution due to the up and down 
movement of the isentropic surface in time (compare with 
Figures 3 and 4a). The mass flux across the upper boundary of 
the lowermost stratosphere was calculated according to 
equation (A8) using UKMO data together with radiatively 
derived mass streamfunction y based on UARS constituent 
measurements [ Rosenlof , 1995]. 
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